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TECHNICAL  NOTE  2^51 

MATHEMATICAL  IMPROVEMENT  OF  METHOD  FOR  COMPUTING  POISSON 
INTEGRALS  INVOLVED  IN  DETERMINATION  OF  VELOCITY 
DISTRIBUTION  ON  AIRFOILS 
By  I.  Fliigge-Lotz 

SUMMARY 


The  Poisson  integral  involved  in  the  determination  of  the  change 
in  velocity  distribution  resulting  from  a  change  in  airfoil  profile  in 
parallel  incompressible  flow  is  solved. 

First,, three  well- developed  numerical  methods  of  evaluating  this 
integral,  all  based  on  the  division  of  the  range  of  integration  into 
small  equal  intervals,  and  the  difficulties  involved  in  each  method, 
are  discussed.  Then  a  new  method,  based  on  the  use  of  unequal  intervals, 
is  developed,  and  compared  with  the  other  methods  by  means  of  several 
examples.  The  new  method  is  found  to  give  good  results  for  both  the 
direct  and  inverse  airfoil  problems  and  is  easily  adaptable  to  rather 
complicated  problems.  It  is  particularly  recommended  for  all  those 
functions  where  steep  slopes  in  small  portions  of  the  region  to  be 
integrated  exist. 


INTRODUCTION 


The  ordinary  thin  airfoil  at  small  angles  of  attack  produces  only 
slight  disturbances  in  the  flow  of  a  parallel  incompressible  fluid. 
Hence,  the  influences  of  camber  and  thickness  upon  the  velocity  distri¬ 
bution  may  be  computed  independently  and  their  effects  superimposed. 

The  effect  of  camber  may  be  represented  by  vortices  distributed  along 
the  chord  line  of  the  airfoil  section;  the  effect  of  the  thickness,  by 
sources  and  sinks  also  along  the  same  chord  line.  The  velocities  pro¬ 
duced  by  these  singularity  distributions  enable  one  to  compute  the 
pressure  distribution  on  the  airfoil  rather  quickly. 

Allen  (reference  1)  has  presented  this  singularity  method  in  a 
form  which  has  proved  to  be  very  practical  for  common  usage.  However, 
in  special  cases  the  unavoidable  evaluation  of  the  Poisson  integral  in 


2 


NACA  TN  2451 


A 


the  course  of  the  computations  has  given  rise  to  numerical  difficulties. 

Such  integrals  are  usually  computed  hy  the  application  of  finite  differ-  / 

ences  using  intervals  of  equal  length.  However,  changes  in  airfoil 
shape,  which  result  in  marked  changes  in  the  function  to  he  integrated 
in  only  small  portions  of  the  range  of  integration,  require  that  extremely 
small  interval  sizes  he  employed  in  this  range,  and,  consequently  over 
the  entire  range  of  integration.  This  leads  to  a  considerahle  amount 
of  computational  work;  hence,  it  appears  reasonable  to  discuss  the 
possibility  of  employing  intervals  of  varying  lengths  for  the  evaluation 
of  the  Poisson  integral. 

This  investigation  was  prompted  hy  the  difficulties  arising  from 
the  problem  of  small  changes  in  the  shape  of  symmetrical  airfoils  at 
the  angle  of  zero  lift.  The  examples  included  in  the  present  report 
are  restricted  to  this  case,  hut  the  results  obtained  are  in  no  way 
specialized  and  may  be  applied  to  all  problems  wherein  the  Poisson 
integral  occurs. 

This  work  was  done  at  Stanford  University  under  the  sponsorship 
and  with  the  financial  assistance  of  the  National  Advisory  Committee 
for  Aeronautics.  ’ 

The  author  wishes  to  express  her  appreciation  to  Mr.  H.  Norman 
Abramson  for  his  intelligent  and  skillful  help  in  the  computational 
work  and  for  his  assistance  in  writing  the  final  report.  The  author 
also  wishes  to  extend  her  thanks  to  Mr.  R.  E.  Dannenberg  and  the 
computing  staff  of  the  J-  by  10-foot  wind-tunnel  section  of  the  Ames 
Aeronautical  Laboratory,  Moffett  Field,  California,  for  preparing  the 
extended  tables  of  the  functions  jno  and  jn0*. 


DISCUSSION  OF  PROBLEM 


The  basic  reference  profile  may  be  given  by  ytr  =  f(x),  and  its 

velocity  distribution  may  be  known  from  an  earlier  computation.  The 
problem  at  hand  is  that  of  determining  the  change  in  the  velocity  dis¬ 
tribution  resulting  from  a  change  in  the  shape  of  the  profile  (indicated 
by  the  dotted  line  in  the  following  fig. ) .  The  difference  of  these  two 
shapes  is  designated  as  Ay^.. 
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Allen  (reference  1,  p.  7)  gives  for  the  change  of  velocity  the  equation 


Av(x0) 

vo 


(1) 


where  VD  is  the  velocity  of  the  basic  parallel  flow.  If,  by 
conformal  mapping  of  the  outside  flow  region,  the  center  line  of  the 
profile  is  transformed  into  a  circle  by  the  relation 


x  =  |  (1  -  cos  9)  (2) 

then  the  profile  is  transformed  into  a  curve  approximating  the  circle 
shown  below. 


I 
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The  change  in  velocity  due  to  a  change  in  form  will  then  he  given  as 


,  __L  f2’  a(Ayt) 


9-9, 


v0  (e°) 


defining 


This  is  the  form  most  often  used  for  computation  purposes  because  the 
inverse  problem  (that  of  computing  the  change  in  shape  due  to  a  change 
in  velocity  distribution)  utilizes  the  analytic  form 


d  (Ay-, 


cot 

V„ 


9  -  9. 


defining 


=  (to.) 

jt+0  V  O/ jr-0 


which  is  strikingly  similar.  The  corresponding  formula  in  the  original 
x.y  coordinate  system  is  given  by 


_  dx  _ 


0  x  - 


k(c  -  x) 


The  evaluation  of  equation  (3)  may  be  accomplished  by  any  one  of 
several  different  methods;  however,  all  of  these  methods  employ  the 
device  of  replacing  the  integral  over  the  range  0  to  2 n  by  a  sum  of 
integrals  over  intervals  of  equal  length  A0.  The  equally  distributed 
points  9n  have  corresponding  values  Xjj  which  are  not  equally 

distributed  (see  following  fig.). 
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e3  01*- 


This  arrangement  is  sometimes  favorable,  and  sometimes  not,  depending 


upon  the  particular  form  of 
equation  (4-3).) 


.  (See  discussion  following 


The  use  of  the  angular  coordinate  9  has  the  advantage  that  the 


functions 


or  &  are  periodic  functions  in  9,  and  this 
vo 


periodicity  facilitates  the  organization  of  the  numerical  computations. 
The  disadvantage  arises  from  the  fact  that  these  functions  are  usually 
given  as  functions  of  x,  and,  since  the  analytic  form  is  not  usually 
known  ^  any  transformations  made  will  lead  to  small  errors .  For  example 
d(Ayt)  Av  * 

if  dx  or  y2,  is  given  at  special  points  which  do  not  correspond 

to  9n  =  nA0,  then  the  computor  must  obtain  the  values  of  these  functions 
for  the  values  9- j,  9^,  and  so  forth  by  interpolation. 


DISCUSSION  OF  SOME  OF  THE  EXISTING  NUMERICAL 


SOLUTIONS  OF  POISSON  INTEGRAL 


The  difficulty  encountered  in  the  solution  of  the  Poisson  integral 

arises  from  the  fact  that  the  term  cot  L°  6o)  or  1  (e  _ 

2  x  -  xD 

tions  (1)  and  (3),  e.g. )  approaches  infinity  when  9  approaches  6  or 
when  x  approaches  xQ.  The  difficulty  is  of  much  less  consequence  when 


the  function  -A_X  or  —  is  given  analytically  than  when  a  numerical 

o 

computation  is  undertaken.  As  a  consequence,  any  simple  integration, 
performed  by  replacing  the  integral  with  a  summation  over  smaller  intervals 
always  requires  that  the  interval  in  which  9Q  or  xQ  is  located  be  given'' 
special  consideration  (see  following  fig.). 
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A  majority  of  the  solutions  currently  in  use  have  been  developed 
to  such  an  extent  that,  for  example,  ^(0O)  13  glven  a  sum  of 

products  of  single  values  of  psv  and  known  factors  AQ;  that  It 


and  known  factors  AQ;  that  is. 


v 
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The  coefficients  An0  depend,  upon  the  particular  method  of 

d  (Ay. 

numerical  integration  which  is  employed.  If,  for  example,  — * — : 


replaced  hy  a  step-curve,  that  is,  assumed  constant  in  every  interval 
(see  fig.  below),  one  set  of  values  of  Ano  would  be  obtained. 


replaced  by  straight-line  segments  (see  fig.  below),  in  which  case  a 
second  set  of  values  of  Ano  would  be  obtained. 


of  segments  of  parabolas,  and  so  forth.  Since  the  accuracy  of  the 

resulting  values  of  depends  upon  both  the  character  of  the  approxi- 

vo 

mate  curve  and  the  size  of  interval  taken,  it  is  apparent  that  the  same 
degree  of  accuracy  might  be  achieved  from  many  different  combinations 


dx 


of  interval  sizes  and  approximations  to  the  function 
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Naiman  (reference  2)  has  used  Simpson's  rule  for  computing  the 
Poisson  integral,  which  corresponds  to  the  replacement  of  the  product 

cot  —  hy  segments  of  parabolas.  The  "critical  interval" 

L  a*  2J  *  v 

(  i.e.,  where  cot  - - >  °o)  was  carefully  treated  by  using  differences 

V  /  2  '  d(Ay.)\ 

of  higher  order  I  including  the  fifth  derivative  of  —  — — - I.  Naiman 

divided  the  period  of  2it  in  20  or  40  intervals  and  calculated  the 
corresponding  sets  of  values  of  AnQ.  Other  workers  at  the  NACA  have 

extended  the  calculation  of  these  values  of  A^  to  80  and  l60  intervals 

(unpublished  information).'1' 

Obviously,  the  time  required  for  computing  increases  with  the 

Vo 

number  of  intervals  taken  because  of  the  increased  number  of  multipli¬ 
cations  to  be  performed.  In  addition,  greater  preparations  for  the 
computing  proceaa  are  necessarily  involyed,  particularly  since  the 

-  ~-r  d(Ayt)  needed  must  usually  be  obtained  by  interpolation. 


values  of 


This  interpolation  has  to  be  done  rather  carefully  as  it  is  often  not 
sufficient  simply  to  take  the  values  of  the  plotted  curve  of 
<l(Ayt) 

-  .  This  curve  should  he  checked  hy  difference  tables  if  the 


values 


are  to  represent  a  smooth  curve 


For  those  functions  of  Ayt  which  may  he  well-represented  hy  a 
Fourier  series,  there  exists  a  simple  method  of  evaluating  the  Poisson 
integral  which  has  apparently  been  overlooked  until  the  present  time. 
This  method  has  the  advantage  of  leading  to  a  computation  which  does 
not  involve  the  derivative  of  Ay^. 


Caiman  has  also  suggested  a  second  method  for  computing  the 
Poisson  integral  (see  reference  3).  In  this  second  method  he  uses 

d(Ay+) 

Fourier  polynomials  to  represent  the  function  — - ^  .  The  computing 


procedure  is  very  simple;  however,  the  results  depend  largely  on  the 

d(Ayt) 

degree  of  approximation  of  — ^ — <-  hy  such  a  polynomial.  Thus,  for 

CLX. 

large  families  of  functions  results  are  good;  however,  cases  are  known 
to  the  author  where  results  were  not  satisfactory  because  regions  with 
steep  gradients  may  not  be  represented  well  enough  by  a  Fourier 
polynomial  of  moderate  order. 
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Equation  (6)  may  be  written  in  a  different  form  (reference  1, 
equation  (43) )  as  follows: 


Av  _  1  rdK)  sin  e 
Vo  Jtjo  dx  cos  0  -  cos 


(8) 


and  this  may  he  rewritten  as 


m  =  12  r  d(Ayt)  IQ  (9) 

Vo  jt  c  Jo  d0  cos  0  -  cos 

Equation  (9)  is  strikingly  similar  to  an  integral  ocurring  in  the 
theory  of  the  lift  distribution  of  a  finite  wing  in  incompressible  flow. 
There  y  the  induced  angle  ai  is  given  by 


a. 


1  f*  d7 _ d0J _ 

2jt  jo  d0f  cos  0 1  -  cos  0 


(10) 


where  7  is  the  local  dimensionless  circulation. 

Multhopp  (reference  k)  has  given  a  solution  for  equation  (10).  He 
divides  the  range  of  integration  into  ^m^  +  1^  intervals  (see  fig. 

below) 


with 


9n  = 


n 


m,  +  1 


yn  =  ?(e n) 


(11) 
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and  computes  at  the  points  0n.  He  assumes  that  y  may  he 

expanded  in  the  form 

7  =X!Cki  sin  V-9 


or 


7 


2 

m  +  1 


S  sin  ^n  sin 

n=l  n=l 


(12) 


The  prime  on  the  summation  symbol  indicates  that  n  =  v  is  to  he 
omitted  from  the  summation  because  that  special  term  has  already  been 
considered  in  the  first  term  of  the  right-hand  side  (i.e.,  bvv7v)* 
Reference  4  presents  tables  for  the  coefficients  bvv  and  bvn  for 
=  7,  15,  and  31.  Applied  to  the  problem  at  hand,  m^  =  31  would 
appear  to  be  rather  small]  therefore  a  table  for  m^  =  63  has  been 

computed  and  is  included  in  the  present  report  (appendix  A).  As  a 
comparison:  For  =  63,  Ad  =  2.8125°;  for  Naiman*s  method  with 

l60  points,  A9  =  2.25°. 

Utilizing  this  method  of  integration  which  was  developed  by  means 
of  Fourier  series,  an  expression  may  be  obtained  for  the  velocity 
distribution  as  follows: 


The  great  advantage  of  this  method  is  that  of  simplicity:  (l)  The 
actual  computational  procedure  is  very  simple  and  (2)  the  derivative 


d(^yt)  is  avoid_edt  The  simplicity  of  computation  is  reflected  in  the 
dx 

fact  that  the  time  required  for  computing  ^  at  one  value  of  dQ 


is  approximately  half  that  required  by  the  method  of  Naiman  when  the 
intervals  have  approximately  the  same  size.  It  should  be  noted. 
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however,  that  the  accuracy  of  the  method  of  Naiman  will  he  greater  than 
that  of  Multhopp  in  those  cases  where  the  differentiation  of  Ayt  hy 
Fourier  expansion  (equation  (12))  does  not  give  good  results. 

A  third  method  of  evaluating  the  Poisson  integral  became  known 
during  the  course  of  the  present  investigation.  In  a  paper  by  Timman 
(reference  5),  the  integral  is  studied  in  the  form 


(15) 


Timman  assumes  that  cr (ijr )  is  not  given  analytically,  but  only  at 
equidistant  points.  An  interpolation  polynomial  (reference  6)  for  b'(\|r) 
is  employed,  and  these  polynomials  replace  the  function  a(\|r)  in  a 
single  interval  by  a  function  of  third  order.  The  polynomial  function 
thus  introduced  has  a  continuous  first  derivative,2  and  it  is  evident 
that  this  continuity  is  essential  for  the  attainment  of  good  results. 

Timman  has  divided  the  period  2ir  into  36  intervals  of  equal 
length  and  established  a  computing  scheme.  The  function  or(iJ/ )  is 
separated  into  its  symmetrical  and  unsymmetrical  parts  so  that 

ffW  =  s  +  d  (16) 

Then 


18  18 

<17> 

where  the  factors  and  PkZ  are  given  in  tabular  form.  In  the 

d(Ay+) 

present  particular  case  — * — in  is  antisymmetrical  (equation  (3))  and 
A  ^ 

~  is  symmetrical  (equation  (4)).  Thus  the  separation  indicated  by 
vo 

equation  (l6)  does  not  require  any  additional  work. 


Timman ’s  method  should  give  good  results  provided  that  a  sufficient 
number  of  intervals  are  taken  -  the  division  of  36  intervals  over  a 
period  of  2jt  (i.e.,  18  intervals  over  the  chord  of  the  profile)  appears 


to  be  insufficient  for  an  accurate 


occurs 


ml 


representation  of  the  function  which 


^The  polynomials  used  in  the  classical  interpolation  formulas  are 
less  smooth  (see  reference  pp.  7  and  10,  figs.  1  and  2). 


12 


NACA  TN  2451 


The  time  required  for  computing  one  point  by  the  method  of  Timman 
is  approximately  the  same  as  for  Naiman's  method  with  the  same  interval 
size. 


Other  methods  of  evaluating  the  Poisson  integral  have  been . 
suggested.  They  will  not  be  discussed  here  as  it  is  the  intention  of 
this  section  to  consider  only  the  most  practical  of  the  known  methods. 
The  three  methods  already  discussed  have  their  own  particular, 
advantages  and  have  been  especially  developed  for  rapid  and  simple 

computation;  however,  all  three  of  these  methods,  when  ^  or  ^  } 

change  rapidly  in  magnitude,  become  cumbersome,  and  require  that  very 
small  intervals  be  taken  over  the  entire  range  of  integration  because 
the  scheme  of  equal  interval  size  is  utilized. 


EVALUATION  OF  POISSON  INTEGRAL  BY  A  METHOD 
EMPLOYING  UNEQUAL  INTERVALS 
Development  of  Method 


As  the  change  in  airfoil  shape,  or  the  change  in  velocity  distri¬ 
bution,  is  given  originally  as  a  function  of  x  it  appears  logical  to 
retain  the  coordinate  x  in  selecting  the  size  of  the  different 
intervals.  Hence,  the  Poisson  integral  may  be  studied  in  the  form 


T 


a(x) 


(18) 


which  corresponds  to  equation  (l).  Conforming  with  its  physical 

meaning  cr(x)  =  ■  -  — —  is  assumed  to  be  a  function  which  is  finite  in 
°  dx  o 

every  point)  of*  its  range  of*  definition. 


Define 

**n  =  *n+i  -  *n 

with 


(19a) 


n  =  0,  1,  2,  3j  •  •  • 


3This  restriction  will  be  dropped  later;  see  discussion  beginning 
with  the  first  paragraph  after  equation  (32). 
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and 


<  x0  <  x  _ 
m  o  ^  m+1 


For  convenience,  there  is  chosen  (see  following  fig. ) 


_  xm  +  xm+x 


— ....  4  t - 1 _ 1 _ 1 _ 1 _ , _ 

•  1  < 

'  1  771 

V  C 

xm 

X 

/  \  . 

xm+l 

O 

The  function  a(x)  is  approximated  by  straight-line  segments 
sketch  in  preceding  section) .  Then,  for  xn  <  x  <  x^j, 

cr(x)  =  cr(xn)  +  n+1^  — ^I^x  - 


=  CTn  + 


an+l  - 


Ax 


n 


“(x  -  *n) 


(19b) 


(19c) 


(see  third 


(20) 


from  which  there  is  obtained 
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Introducing  j  into  equation  (21),  there  results 


™  ■- j-X!( 


an+l  -  °n)  + 


cr  -  a 
n+1  n 


^(Xo  -  *n) 


i  5>Jn°  +  2(ffn+l  "  CTJ  f1  + 


xo  -  *n  . 
Axn  Jno 


Or,  defining 


i  +  *n  , 

Axn  Jno 


there  results,  finally. 


(xo)  = 


anJno  +  /*  |  (an+l  “  an)  <)no* 


Since 
written  as 


xn+l  =  ^  +  &Xn>  the  functions  jno  and  jnQ*  may  be 


l  *  =  i  +  xq  -  j 

'no  ay  Jno 


Jno  =  log^l1  + 


=  log, 


xn  “  xo 


<-i  + 


xo  “  xn 


for  ^  >  x( 


for  xn  +  Axn  >  xQ  >  ^ 


=  1°ge\v1  +  x^~ -  xj  for  xo  >  ^  + 


and  this  form  shows  that  jnn  and  j  *  are  functions  of  Xl1  ~  X° 


(27a) 


only, 
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~  XQ  _ v.  ± 


For  xQ  -  xn  =  |  Axn 


For  very  large 


Jnn  =  0 


6, 


Jno  I  -  i  +  • 


For  very  large  negative 


xn  "  xo  .  xn  xo 
-  with  - 


-  ft*. 


_1 _ 1_ 

,*  /.*2 

I  26 


1  *  _>  _  _1 _ L_  + 

Jno  '  *  p 

21  36* 


(2Tb) 


(27c) 


(27d) 


These  functions  are  given  in  figure  1  and  in  table  I. 

It  is  seen  that  -high  absolute  values  of  jno  and  jno  occur 

x  -  xQ  _  4 

near  those  values  of  which  characterize  the  critical  interval* 

■^n 

Figure  1  gives  an  idea  of  the  characteristic  qualities  of  jno 

-  .  %  "  XQ  _  - - lUrv  norurtoomrf  of  1  1  <5  nnt. 


and  jn0  as  functions  of 


•  however,  the  representation  is  not 


sufficient  for  picking  out  values  for  a  computation.  Table  I  gives  the 

X  - 

■X  _  i  ^  XI  U  1.  r\  r-  nru  •  _  ^  m  -T 


values  of  jnQ  and  jno  for  -1+9*5  < 


<  1+9*5*  This  table  might 


Xq  =  *  5°±1  (equation  (19c))  the  critical  interval  i 


<  0.5- 


given  by  -0*5  < 
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be  used  for  rough  computation  and  for  getting  acquainted  with  the 
method.  In  general,  it  is  advisable  to  use  those  tables  which  are 
given  in  appendix  B. 

It  will  prove  of  benefit  to  investigate  the  exactness  of  that 
portion  of  the  integral  which  contains  the  singularity.  Recalling 
that  the  function  u(x)  was  replaced  by  a  straight  line  in  every 
interval  (equation  (20)),  there  is  obtained: 


ct(x) 


X  -  X, 


=  -  ~(°m+l  -  am) 


(28) 


if  Xjh  <  xQ  <  xm+1<  Now,  let  an  expansion  of  the  function  cr(x)  in 
the  critical  interval  around  xQ  be  assumed  as  follows: 

cr(x)  =  c(x0)  +  a»(x0)(x  -  x0)  +  ^£^(x  -  xQ)2  + 


g,,,(xo) 


3  .  glV(xo) 


31 


7^(x-xo)  +  0Jix  ~  xo) 


4*  .  .  • 


(29) 


Then, 


dx  =  -  - 


0.K)te  +  ^f(f)3  +  .  . 


It 


22(am+l  ”  °m)  ”  ^6  (CTm+2  "  ffni-l) 


36 


‘  n  \^(>+1  “  °m)  ' 

(°m+2  “  am+l)  +  (am  “ 


_1. 

36 


(30) 


Comparison  of  formulas  (30)  and  (28)  shows  that  the  error  in  the 
critical  interval  is  approximately  given  by 
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The  error  of  evaluating  the  whole  integral  by  finite  differences  may 
he  estimated  by  using  two  different  interval  distributions  and  comparing 
the  results  for  a  given  xD. 

However,  in  addition  to  that  error  of  the  result  produced  by 
replacing  the  Poisson  integral  by  a  sum  there  exists  another  error. 

This  sum  cannot  be  computed  exactly,  but  has  a  certain  error  depending 
on  the  accuracy  of  the  given  data  for  crn(x)  and  the  tabulated  values 

*  H&t) 

of  jno  and  jn0  .  As  the  function  a(x)  =  — - —  usually  has  an 

-1  ^ 

error  of  ^  =  1  X  10  J  it  has  proved  amply  satisfactory  to  give  jno 
and  jn0*  to  four  decimal  places,  the  error  being  less  than 
e2  =  5  X  lO-^.  The  error  of 


is  smaller  than  its  upper  limit  given  by 


2lJ. 


2  M +  2 


This  formula  shows  that  the  influence  of  €]_  is  stronger  than  the 
influence  of  e2  as  long  as  2,  1  jcrn  |  +  '|°n+l  "  °n|  is  smaller 

than  1  -  as  it  is  in  our  later  examples  -  and  the  sums 

and  ^Z!|jno*|  are  always  larger  than  1.  An  increase  of  subdivisions 

makes  the  sums  in  the  upper  limit  of  the  error  (32)  grow,  thus  requiring 
a  higher  accuracy,  especially  in  crn  and  perhaps  also  in  the  values 

of  jno  and  jno*. 


In  establishing  the  solution  of  equation  (18)  it  was  assumed 
that  cr(x)  is  finite  throughout  its  range  of  definition.  If  it  is 
desired  to  compute  the  change  in  shape  due  to  a  proposed  change  of 
velocity  distribution,  this  restriction  must  be  eliminated,  as  will 
be  recognized  immediately. 
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% 


* 


Equation  (5)  may  be  written  in  the  form 


dfot) 

dx 


Av/V0  dx 

\/x(c  -  x)  x  -  xQ 


(5’) 


Omitting  the  factor  \J  xQ\c  -  xQ^ ,  which  does  not  affect  the 

process,  the  integral  may  be  reduced  to  the  form  of  equation 
defining 


integration 
(18)  by 


Av/Vp 
\/x(c  -  x) 


Ol(x) 


(33) 


However,  cr^x) 
and  i  0, 


will  be  infinite  at  x  =  0  and  x  =  c  if  ^  0 

Vo/o 

therefore,  a  special  consideration  of  the  neighborhood 


of  x  =  0  and  x  =  c  is  required.  This  is  done  by  splitting  the 
integral  into  the  following  three  parts: 


% 


with  and  being  small  compared  with  c.  The  integral 


may  be  treated  as  was  explained  formerly  for 


dx 


-  x 


o 
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(see  equation  (l8))  because  ^(x)  is  finite  for  <  x  <  c  - 

For  the  first  and  third  integrals,  however,  a  new  integration  formula 
must  "be  developed.  By  introducing 


\1  =  c  -  x 


and  °]_(x)  ~  cri^J'(x)  =  a±  (l1) 


there  is  obtained 


0- I  (x)' 


dx 


X  -  X, 


2  */  \  dM- 

crx  (^) - 

0  ^  -  ^o 


(35) 


c-e, 


Hence,  the  method  used  for  the  first  integral  will  also  apply  to  the 
third.  In  most  cases  )  will  he  zero  and  there  will  he  no  need 

Vole 

for  a  special  evaluation  in  the  neighborhood  of  x  =  c. 


The  integral 


Av/V0 _ dx 

\/x(c  -  x)  x  ■  xo 


(36) 


* 


w 


will  have  an  important  influence  on  the  result  of  equation  (34)  only 
if  xQ  is  near  to  gq_.  First  the  general  formula  will  he  given  and 

then  a  simplification  will  he  discussed  for  xQ  » 

The  integral  (36)  will  he  solved  assuming  that 

^  =  ao  +  al(f)  +  a2(f)  f0r  0  <  X  <  €1  (3' 

Only  the  final  formula  of  this  procedure  is  given  here;  the 
details  of  the  solution  will  he  found  in  appendix  C. 
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* 


1* 


with  M0  given  in  figure  2,  and 


with 


with 


(39) 


The  coefficients  aQ,  a^,  and  a^  may  he  determined  first,  as 
they  do  not  depend  upon  the  particular  value  of  xQ,  and  then  F^x0^ 
may  he  computed.  The  term  depending  on  a^  and  a^  will  exert  an 
influence  only  for  small  values  of  xQy/c.  After  a  brief  training  the 
computor  should  he  able  to  decide  rather  accurately  when  the  formula 


is  sufficient  and  when  the  more  exact  expression  (equation  (38))  is 
required  (also  see  fig.  2). 


V 
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Organization  of  Computational  Procedure  for  Unequal  Intervals; 

Transition  from  One  Size  of  Interval  to  Another 

A  thorough  understanding  of  the  method  is  test  achieved  by 
following  through  a  rather  simple  example;  in  addition  various  short 
cuts  to  the  method  will  he  demonstrated. 

Assume  a  function  cr(x)  of  the  type  shown  in  the  following 
figure . 


It  appears  reasonable  to  take  rather  small  intervals  for  small  values 
of  x  because  of  the  form  of  the  curve  cr(x);  therefore,  the  following 
arrangement  of  interval  sizes  is  arbitrarily  selected: 


Ax  =  0.002  for  0  <  x  <  0.030 

Ax  =  0.006  for  0.030  <  x  <  O.O96 

Compute  t (0.009)  with  the  help  of  equation  (26).  Note  that  the 
critical  interval  extends  from  0.008  to  0.010.  Table  II (a)  gives 
the  values  of  x/c,  an,  an+1  -  crn,  Jn£),  and  jno*  for  the  range 

with  Ax  =  0.002.  At  x  =  0.030  the  interval  changes  to  Ax  =  0.006 
and  the  same  functions  are  given  for  the  range  with  this  size  of 
interval  in  table  II (b).  Naturally  the  range  above  the  broken  line 
in  table  II (b)  is  not  utilized  in  the  computation  since  this  portion 
has  been  considered  in  table  II (a). 
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Note  that 


progresses  in  table  II (b)  in  the  same  manner 


as  in  table  II (a);  this  is  due  to  .the  special  choice  of  Ax,  If 

Ax  =  0.006  were  used  starting  with  x  =  0  the  critical  interval  for 
xQ  =  0.009  would  extend  from  0.006  to  0.012,  Hence,  for  Ax  =  0.006, 

jno  =  0  and  jno  =1  is  to  be  found  at  x/c  =  0,006. 


For  rapid  computation  it  is  best  to  have  jnQ  and  <jno*  as 

xn  "  xn 

functions  of  ^  ■■  on  a  paper  strip  and  to  place  this  strip  adjacent 

x  —  x 

to  the  columns  headed  by  crn  and  an+1  -  an.  If  ----- — —  progresses 

ZAX 

as  indicated  in  table  I,  the  correct  location  of  jn0  =  0  and  jno*  =  1 

at  the  beginning  of  the  critical  interval  fixes  the  placement  of  the 
strip. 


In  the  example  just  treated,  the  transition  from  one  size  of 
interval  to  another  is  very  easy  because  xQ  lies  at  the  midpoint  of 
an  interval  of  the  size  0.00 6  as  well  as  of  the  size  0.002,  if  starting 
with  x  =  0. 

If  Ax  had  been  chosen  0.004,  such  a  desirable  arrangement  would 
not  have  resulted  because  xQ  =  0.009  would  not  be  located  at  the  mid¬ 
point  of  an  interval  of  this  size  (starting  with  such  intervals 
at  x  =  0 ) . 


As  a  second  example  compute  the  value  of  t  at  xG  =  0.015. 

*n  ~  xo  *n  "  xo 

Again,  - = —  and  - — —  will  progress  as  in  table  I.  The  values 

Ax  Ax 

jno  “  0  an<3-  4io*  =  1  will  be  placed  opposite  x/c  =  0.014  for  the 

region  with  Ax  =  0.002  and  opposite  x/c  =  0.012  for  the  region 

with  Ax  =  0.006.  As  long  as  Ax  =  3Ax,  Ax  =  3Ax,  and  so  forth  and 
if  xQ  is  chosen  so  as  to  be  at  the  midpoint  of  the  largest  size  of 
interval,  the  computation  may  be  accomplished  by  shifting  the  strip 
with  jn0  and  jnQ*  corresponding  to  table  I. 


But  suppose  that  the  interval  sizes  are  so  arranged  and  it  is 
desired  to  compute  a  point  where  xQ  does  not  lie  at  the  midpoint  of 

the  largest  size  of  interval;  for  example,  xQ  =  0.013.  The  value 

xQ  as  0.013  lies  at  the  midpoint  of  an  interval  with  Ax  =  0.002;  hence, 


24 


NACA  TN  2451 


for  the  range  0  <  x  <  0.030,  jno  and  jno*  may  be  taken  directly 

from  table  I.  However,  at  x/c  =  0.030,  intervals  of  the  size 
Ax  =  0.006  commence  and  there  is  obtained 


*n  -  Xo 


Sc 


0.030  -  0.013 

0.006 


2.833 


The  value  of  — - —  progresses  by  1,  that  is,  2.833>  3.833* 

Ax 

4.833,  ....  Thus  the  functions  jno  and  jno  are  needed  for 

values  of  — - —  which  are  not  given  in  table  I.  One  might  think  of 

Ax 

taking  them  out  of  an  enlarged  diagram  (see  fig.  l);  however,  it  is 
much  more  convenient  to  take  them  out  of  an  extended  table,  which  is 
conveniently  arranged  for  "advancing  by  1. "  Such  tables  are  given  in 
appendix  B. 

The  example  presented  by  the  figure  at  the  beginning  of  this 
section  suggested  starting  at  x  =  0  with  the  smallest  intervals. 
However,  other  examples  may  suggest  another  distribution  of  intervals. 
The  smallest  size  of  intervals  may  lie  at  any  part  of  0  <  x  <  c. 

There  are  no  restrictions  in  the  arrangement  of  intervals.  (See,  e.g. , 
discussion  following  equation  (43).) 


¥ 


Accuracy  of  Method,  Examined  by  Means  of 
an  Analytical  Example 

The  accuracy  of  the  result  depends  directly  upon  the  size  of  the 
interval  taken  and  the  reliability  of  the  data  comprising  the  func¬ 
tion  u(x).  Because  the  function  a(x)  will  be  replaced  by  a  broken  line, 
a  glance  at  the  curve  will  quickly  suggest  an  arrangement  of  intervals. 

In  addition,  the  error  in  the  critical  interval  may  be  used  as  a  first 
test  of  the  choice  of  intervals. 

As  a  test  of  the  quality  of  this  new  method,  involving  unequal 

d(Ay.) 

intervals,  a  function  cr(x)  =  — -  has  been  treated  which  allows 

dx 

the  analytical  computation  of  t (x)  =  ^  . 


t 
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The  function  cr(x)  is  given  analytically  as 


0  Is  x  <  2A 


d(Ayt)  _ 


dx 


=  Bx(2^  -  x) 


2A  =  x  =  cx 


=  -D(oi  -  x)(x  - 


2A) 


C1<X<C 
1  dx 


The  following  arbitrary  values  have  been  selected: 


c±  =  0.35 


c  =  1.0 


1  d(Ayt) 

D  =  Some  multiple  of  B  so  that  /  — - - -  dx  =  0 


0 


dx 


dfot) 


The  functions  Ay.  and  — - — ^  are  given  in  figures  3(a)  and  3(h), 
x  dx 

respectively. 
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The  analytical  computation  of  ^  for  figure  3(h)  is  given  in 

vo 

figures  4(a)  and  4(h).  The  arrangement  of  the  unequal  division  for  the 

numerical  computation  of  is  indicated  in  figure  4(a). ^ 

*  o 

^ It  was  desirable  to  obtain  the  value  of  ^  at  xQ  =  0;  hence , 

vo 

the  first  interval  has  been  placed  so  that  -0.001  <  x0  <  0.001  and  the 

^(^t) 

function  ^  =  0  for  -0.001  <  x  <  0.  Since  the  function  -  - 

dx  # 

is  replaced  in  every  interval  by  a  straight  line,  the  error  might  be 
expected  to  be  large.  However,  g  =  0  at  Xq  =  0  will  aid  in 
preventing  the  error  from  being  too  large. 


A 

g! 

y  / 

s  / 

\ 

\ 

\ 

_  Ax 

0 

o 

A  more  exact  solution  would  be  obtained  by  putting 


g  =  0  -f<*<0 

g  =  Sx  0<X<^ 
Ax  ^ 


1  1 


«J0  )*-*o  *  " 

For  the  interval  -  —  <  x  equation  (26)  would  yield 


PAx/2 
-  Ax/2 


X  -  X. 


-  ^  -  0  )  Jno  +  0  X  jno 


and  no  error  is  introduced.  For  xQ  ^  0  there  is  a  very  small  error 
which  may  be  avoided  by  respecting  the  change  of  size  of  the  interval 
near  x  =  0. 
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* 


Also  given  in  figure  4(a)  are  points  of  the  ^  curve  determined  by 

*o 

the  method  of  unequal  intervals.  Figure  4(b)  presents  the  same 
information  plotted  to  a  larger  scale. 

For  comparative  purposes  the  same  problem  has  been  treated  by  the 
three  methods  of  computation  discussed  earlier,  namely,  those  of  Naiman, 
Multhopp,  and  Timman.  Figures  5(a)  and  5(b)  show  the  results  obtained 
by  the  method  of  Naiman;  obviously,  the  40-point  solution  does  not  use 

a  sufficiently  accurate  representation  of  the  )  curve,  while  the 

dx 

80-  and  160-point  solutions  are  quite  good,  with  the  exception  of  the 
maximum  and  minimum  points  of  the  ^  curve.  In  order  to  obtain  a 

value  at  approximately  |  =  O.O36  a  solution  involving  320  points  would 

be  required.  In  this  respect  the  method  of  unequal  intervals  is  more 
adaptable  to  special  conditions  without  involving  much  new  work  than  is 
the  method  of  Naiman. 


The  results  obtained  by  Multhopp’s  method  are  given  in  figures  6(a) 
and  6(b).  The  31-point  solution  (in  Multhopp 's  somewhat  odd  manner  of 
designation)  corresponds  to  A 9  =  5.625°;  the  63-point  solution,  to 
A  B  =  2.8125°.  The  computation  is  very  simple  and  the  results  of  the 
method  with  63  points  are  comparable  with  that  of  Naiman  with  80  points, 
with  the  exception  of  those  near  the  region  0  <  x  <  0. 01  (this  is 

shown  most  clearly  in  fig.  6(b)).  The  very  steep  peak  of  — ^"*^0 

cLx 

at  x/c  =  0.02  requires  rather  high  harmonics  for  the  representation 
of  Ay^;  consequently,  good  accuracy  in  the  region  near  the  origin  may 
not  be  expected.  This  is  substantiated  by  the  fact  that  for  the  63-point 
method  the  highest  effective  harmonic  would  have  three  waves  in  the 
region  0  <  x  <  0.04;  obviously  a  sufficient  degree  of  accuracy  in  the 
differentiation  process  cannot  be  obtained. 


As  mentioned  earlier,  Timman 's  method  might  be  expected  to  give 
good  results  if  the  size  of  interval  is  properly  chosen.  Inasmuch  as 

only  a  table  for  A 9  =  =  10°  was  available,  the  result  of  the 

computation  for  ^jr-  cannot  be  expected  to  be  good,  as  is  evidenced  by 

observing  figure  7.  The  result  obtained  is  comparable  with  that  of 
Multhopp's  15-point  and  Naiman's  40-point  solutions. 


An  excellent  method  of  examining  the  accuracy  of  these  methods 
still  further  is  simply  that  of  solving  the  inverse  problem.  From  the 

—  just  discussed,  values  for  — (f^t) 

»o  dx 


curves  of 


have  been  computed 
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and  are  presented  in  figure  8.  The  method  of  unequal  intervals  gives 
good  results,  indicating  that  the  arrangement  of  intervals  chosen  was 
as  good  for  the  inverse  problem  as  for  the  direct  problem.  It  is 
apparent  that  Naiman's  method  requires  even  smaller  divisions  than 
160  points  in  order  to  avoid  inaccuracies  near  the  point  x/c  =  0.04. 

The  reader  may  wonder  that  the  inverse  problem  is  not  given  by 
Multhopp ' s  method.  It  must  be  recalled  that  Multhopp's  method  of 
solving  the  direct  problem  does  not  involve  the  differentiation  of  Ay^; 

that  is,  it  is  particularly  fit  for  this  problem  and  presents,  on  the 
other  hand,  no  analogy  for  the  inverse  problem: 


d(Ayt)  _  ±  f2*  Av  6  -  0o  d0  =  _  1  f*  AI - 3J±1° - d0 

dx  2n  Jo  VQ  2  « Jo  V0  cos  9  -  cos  90 

Because  more  extended  tables  for  Timman»s  method  are  not  available, 
and  the  results  obtained  from  the  36-point  method  for  which  tables  exist 
are  very  poor,  no  further  examples  of  the  application  of  this  method 
will  be  given. 


COMPARISONS  OF  METHODS  OF  NAIMAN  AND  MULTHOPP  WITH  METHOD  EMPLOYING 
UNEQUAL  INTERVALS  BASED  ON  ACTUAL  EXAMPLES  OF  CHANGES  IN 

AIRFOIL  SHAPE 


The  method  of  unequal  intervals  has  shown  good  qualities  when 
applied  to  a  problem  where  the  function  cr(x)  is  known  analytically. 
However,  as  mentioned  earlier,  this  function  is  not  usually  known  in 
analytic  form.  This  section,  therefore,  will  compare  the  three 
principle  methods,  those  of  Naiman,  Multhopp,  and  unequal  intervals, 
on  the  basis  of  actual  design  problems,  solving  the  direct  problem 

for  —  and  using  these  results  to  solve  the  inverse  problem  (excluding 
Vo 

Multhopp  for  the  inverse  problem). 


Figure  9(a)  shows  the  Ayt  relations  for  examples  I 
d(Ay  ) 

figure  9(b),  the  Xjj —  relations.  Note  that  the  slope 

for  example  II  is  more  than  twice  that  of  example  I  near 


and  II  and 
^t)' 


of  1  dx  J 
x/c  =  0. 
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The  direct  problem  for  example  I  hy  Naimanfs  method  is  given  in 
figure  10.  The  l60-point  solution  does  not  show  any  appreciable 

deviation  from  the  80-point  solutions  at  the  region  of  (  ; 

\  ^Vo/max 

however,  near  the  origin,  at  /— M  ,  the  influence  of  the  smaller-sized 

\Vo/min 

intervals  (8O:A0  =  4.5°;  l6O:A0  =  2.25°)  is  quite  pronounced. 

The  solution  by  Multhoppfs  method  is  given  in  figure  11;^  31  points 
around  the  half  circle  are  not  sufficient  for  a  solution  comparable  with 
Naimanfs  80-point  solution,  and  even  a  solution  based  on  63  points  does 
not  offer  much  improvement.  The  results  are  poor,  as  might  be  expected, 
in  the  region  very  near  the  origin  (see  preceding  section). 

Figure  12  presents  the  results  obtained  by  the  method  of  unequal 
intervals,  compared  with  results  obtained  by  Naimanfs  80-  and  160-point 
solutions.  The  method  of  unequal  intervals  gives  results  corresponding 
to  those  established  by  Naiman's  l60-point  solution.  The  subdivision 
used  is  shown  in  the  figure. 


As  before,  the  inverse  problem  was  solved,  and  is  given  in  fig¬ 
ure  13.  In  each  case  the  computed  curve  of  was  the  one  used  in 


obtaining  the  values  for  the 


d(Ayt) 

dx 


curve.  Both  methods  give  good 


results,  thus  proving  that  the  chosen  number  of  divisions  was  sufficient 
in  Naimanfs  method  and  in  the  method  employing  unequal  intervals. 


The  value  of 


dfot) 

dx 


computed  at  x/c 


0.171  is  of  some  interest. 


This  point  was  computed  by  the  method  of  unequal  intervals  in  two 
different  ways:  First,  the  arrangement  of  intervals  shown  in  figure  13 
was  utilized  to  compute  the  lower  point.  Then  a  new  arrangement  of 
intervals  (Ax  =  0.018  for  0  <  x  <  O.36)  was  set  up  and  the  same  point 
computed.  The  idea  was  to  determine  the  inaccuracies  that  would  result. 
One  might  predict  that,  since  the  point  x/c  =  0.171  lies  at  a  con¬ 
siderable  distance  from  the  region  of  rapid  changes  in  errors  of 

Vo 

only  small  magnitude  would  be  introduced;  this  is  fairly  well  sub¬ 
stantiated  by  the  results  shown  in  the  figure  because  the  error  thus 
introduced  is  approximately  that  of  the  deviation  of  Naiman's  160-point 
solution. 


Recall  that  this  method  does  not  involve  the  differentiation 
of  Ayt. 
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Now,  turning  our  attention  to  example  II  ^  which,  it  will  be 

recalled,  has  a  slope  of  — ^  of  approximately  twice  that  of 

dx 

example  I ,  the  results  given  in  figures  14  to  17  are  obtained. 


For  the  direct  problem  Naiman,s  method  of  l60  points  and  the 
method  employing  unequal  intervals  give  results  which  are  in  good 
agreement.  For  the  inverse  problem  (fig.  17)  it  is  apparent  that  the 


method  using  unequal  intervals  is  superior  (  see  the  deviation  at 

d(4yt) 


dx 


-1  max 


given  by  the  method  of  Naiman^.  Multhopp's  method  gives 


a  rather  good  result  (fig.  15),  which  may  be  attained  when  the  Fourier 
representation  of  the  Ay  curve  is  adequate. 


The  two  examples  thus  far  presented  are  favorable  for  Kaiman's 
method  because  the  steep  slopes  of  cr(x)  occur  near  x  =  0  wh<£re  the 
points  Naiman  uses  are  close  together.  However,  going  to  still  steeper 
slopes  near  x  =  0  would  require  a  rapidly  increasing  number  of  points. 
The  new  method  offers  another  possibility  here.  Assume  that  in  that 
critical  region  x^.  <  x  <  (xk  ^y  0)  d(x)  ^sy  represented 

by  <r(x)  =  E  anxn.  Then  the  integral 


may  be  split  into  three 


integrals 


The  first  and  third  of  these  integrals  may  be  solved  in  the  usual 
manner  using  the  functions  jno  and  second  integral  will 

be  solved  analytically. 
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This  simple  form,  due  to  the  use  of  the  coordinate  x  in  the 
Poisson  integral,  allows  a  rapid  integration,  because  the  integral 


can  he  solved  by  recurrence  as  follows: 


x  n  -  x  n 

=  JS±- - fc-  +  x0kn_!  0  for  n  £  1 

'  r-i  ’ 


(42) 


n 


with 


* 


/ 


(43) 


Thus  even  very  steep  slopes  cause  no  difficulties. 

As  already  mentioned,  examples  I  and  II  correspond  well  to  the 
qualities  demanded  by  Naiman's  method  insofar  as  the  rather  steep 
slopes  occur  in  those  portions  where  the  points  0n  are  close  together. 
If  those  steep  slopes  should  occur  in  other  portions  of  the  chord,  how¬ 
ever,  a  very  great  number  of  points  in  the  Naiman  method  would  be  needed 
in  order  to  represent  cr(x)  adequately,  and  to  get  reliable  results. 

In  such  a  case  the  method  using  unequal  intervals  shows  its  advantage 
by  allowing  a  free  subdivision  of  the  chord. 


A  third  example  will  serve  to  illustrate  this.  Figure  18  shows  a 


function  cj(x) 


=  d(Ayt) 


dx 


The  essential  values  of  the  function  lie  in 


a  part  of  the  chord  where  even  Naiman's  method  with  160  points  is  not 
sufficient  to  represent  the  function  accurately.  This  is  forcibly 

shown  by  the  two  curves  of  If  the  function  o(x)  is  modified 

vo 

(dotted  line)  so  as  to  eliminate  the  high  peak,  then  the  curve  by 

vo 

unequal  intervals  can  be  made  to  agree  with  the  original  ^  by  Naiman's 

vo 

160-point  solution,  thus  definitely  proving  that,  in  this  example, 
Naiman's  method  with  160  points  Is  insufficient. 


Table  III  indicates  the  computation  for  the  point  xQ  =  0.065 
by  unequal  intervals. 
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CONCLUDING  DISCUSSION 


The  new  method  of  evaluating  the  Poisson  integral  developed  herein 
is  to  be  recommended  for  all  those  functions  ct(x)j  where  steep  slopes 
in  small  portions  of  the  region  to  be  integrated  exist.  In  these . 
portions  a  very  small  size  of  interval  may  be  chosen  without  requiring 
that  this  same  size  of  interval  be  used  throughout  the  region  of 
integration.  In  this  manner,  the  work  required  for  computation  may  be 
maintained  at  a  reasonable  level  even  for  the  most  complicated  problems 

The  analytical  treatment  of  special  parts  of  the  integral  is 
possible  (evaluating  the  remainder  by  the  new  method;  see  preceding 
section).  In  those  problems  where  a  transition  to  very  small  intervals 
in  part  of  the  integration  range  would  require  the  determination  of  a 
great  many  values  of  cr^,  this  idea  might  be  used  to  advantage. 

It  should  be  noted  that  the  smoothness  of  the  function  a(x)  and 
its  accurate  representation  by  single  points  is  essential  for  good 
results.  If,  for  example,  single  points  an  are  simply  taken  from  a 
curve  for  x^  very  close  to  one  another  it  may  be  compulsory  to  check 
these  values  by  a  table  of  differences. 
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APPENDIX  B 

VALUES  OF  jno  AMD  Jn(*  AS  FUNCTIONS  OF  Xn~X° 

Ax 

The  values  of  the  functions  jno  and  jno*  are  presented  as  indi¬ 
cated  in  the  following  table.  The  values  are  tabulated  in  a  form  selected 
to  minimize  the  necessity  for  interpolation  except  for  the  region  contain¬ 
ing  the  singularities  of  the  functions  jno  and  j  *.  For  ease  in 
computation,  tables  B— I  to  B— VIII,  inclusive,  are  arranged  so  that  the 

vertical  increment  of  is  unity.  Table  B-IX  gives  additional  values 

for  the  region  containing  the  singularities  of  the  functions  i_  and 
a  *  °  no 

Jno  * 


TABLE  NO. 

RANGE  OF  Xn  x° 

Ax 

INCREMENT  OF  Xn“X° 
Ax: 

B— I 

— 189  to  —90 

1.0 

B— II 

-89.5  to  -40.0 

.5 

B— III 

-39.9  to  -20.0 

.1 

B— IV 

-19.99  to  0 

.01 

B-V 

0  to  19.99 

.01 

B— VI 

20.0  to  39.9 

.1 

B— VII 

4o.o  to  89.5 

.5 

B— VIII 

90  to  I89 

1.0 

B-IX 

-1.000  to  0.000 

0.001 

VALUES  OF  AND  USED  IN  EVALUATING  EQUATION  (26) 
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HrlHHrlHrlHHHHHHrlOJWOJfllWCy 
OOOOOOOOOOOOOOOOOOOO 

1  1  1  1  1  I  1  1  1  1  1  1  1  1  1  1  1  II 

,9  9  onco  t-co  0H4  f-ojcovo^-  iTvco  m 

vovo  r— co  os  oso  h  oj  m  A'S  on  0  cm  ^  vo  co  h 
oj  oj  oj  OJOJOjmmmmmmmm.d'.4--4-_4--4-  A 
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APPENDIX  C 


DETAILS  OF  SOLUTION  OF  INTEGRAL  (36) 


F1  = 


x(c  -  x)  X  -  X. 


1  /  A V_1 

c  Jo  V„  \[ 


]_  .  S  X  „ 

c  c  c 


Introduce  £  =  —  and  find 
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^  Av  _1_  (1  +  1  |  +  i  |2  + 
V-  \/T  \  2  8 


F  _  1  /  Av  1 

1  ~  c7o  VQ 


dl 

I  -  I, 


With  the  expansion  (see  equation  (37)) 


\  -  a0  +  ^al  +  |  =  +  ^a2  +  g  ai  +  8  S°)  ' 


Hence, 


■jf  *  2 

aQ  +  a1  I  +  a2  |  + 


i  r^c  ao + ai 


*t  ^  a  *t  2 


an  +  a,  |  +  aP  | 


nTT  (i  -  lo) 


P€llc 

1  '  dl 

=  ‘  a°J 0  vr(t-  e0) 


*  rei/' 


^Mo) 


■l/C  E  2 


I g  dl 

\l7(t  -  s0) 


56 


NACA  TN  2451 


8 


NACA  TN  2451 


57 


With  these  expressions  the  integral  is  as  follows 


=  j(aoLo 


*T 

+  a-^  +  a2  ^2 


°  Mq  +  M0  +  +  a2*  (\fio 


2^#  +  f 


+  ao  I  \/50  M0  + 


3\ 


Mo(— +  a 

\/Io 


iVw\/c 


VtK  +  W)  +  f 


(C12) 


Av 

The  coefficients  a0,  a-,,  and  ap  of  the  expansion  of  —  are  given 

-L  c.  V  o 

by 


/Av\ 

a°  =  H 

\  *  0/ 


a,  =  — 
1  2e. 


x=0 


-3 


rAv' 
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X=0  X=€ 


x=2e 


a2 - 2 

20 


/Av\  _  ?/Av\  ,  /  Av\ 

Vo/X=0  Wx=El  V  0/ x=2g. 


(C13) 
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TABLE  I.-  VALUES  OF  jn0  AMD  J 


Jno 


FOR  -49.5  < 


Ax 


*n  -  *o 

Ax 

^no 

.  * 

^no 

*11  -  *0 

Ax 

Jno 

-49.5 

-0.0204 

-0.0102  • 

0.5 

1.0986 

-48.5 

-.0208 

-.0101+ 

1.5 

.5108 

-47.5 

-.0213 

-.0107 

2.5 

.3365 

-46.5 

-.0217 

-.0109 

3.5 

.2513 

-45.5 

- . 0222 

-.0112 

4.5 

.2007 

-44.5 

-.0227 

-.0114 

5.5 

.1671 

-43.5 

-.0233 

-.0117 

6.5 

.1431 

-42.5 

-.0238 

-.0120 

7.5 

.1252 

-41.5 

-.0244 

-.0122 

8.5 

.1112 

-40.5 

-.0250 

-.0125 

9.5 

.1001 

-39.5 

-.0256 

-.0129 

10.5 

.0910 

-38.5 

-.0263 

-.0132 

11.5 

.0834 

>37.5 

-.0270 

-.0136 

12.5 

.0770 

-36.5 

-.0278 

-.0139 

13.5 

.0715 

>35.5 

-.0286 

-.0143 

14.5 

.0667 

-34.5 

-.0294 

-.0148 

15.5 

.0625 

-33.5 

-.0303 

-.0153 

16.5 

.0588 

>32.5 

-.0313 

-.0157 

17.5 

.0556 

-31.5 

-.0323 

-.0162 

18.5 

.0526 

-30.5 

-.0333 

-.0167 

19.5 

.0500 

-29.5 

-.0345 

-.0173 

20.5 

.0476 

-28.5 

-.0357 

-.0180 

21.5 

.0455 

-27.5 

-.0370 

-.0186 

22.5 

.0435 

-26.5 

-.0385 

-.0193 

23.5 

.0417 

-25.5 

-.0400 

-.0201 

24.5 

.0400 

-24.5 

-.01+17 

-.0210 

25.5 

.0385 

-23.5 

-.01+35 

-.0219 

26.5 

.0370 

-22.5 

-.01+55 

-.0229 

27.5 

.0357 

-21.5 

-.0476 

-.0240 

28.5 

.0345 

-20.5 

-.0500 

-.0252 

29.5 

.0333 

-19.5 

-.0526 

-.0266 

30.5 

•  0323 

-18.5 

-.0556 

-.0280 

31.5 

.0313 

-17.5 

-.0588 

-.0297 

32.5 

.0303 

-16.5 

-.0625 

-.0316 

33.5 

.0294 

-15.5 

-.0667 

-.0337 

34.5 

.0286 

-14.5 

-.0715 

-.0362 

35.5 

.0278 

-13.5 

-.0770 

-.0390 

36.5 

.0270 

-12.5 

-.0834 

-.0423 

37.5 

.0263 

-11.5 

-.0910 

-.0462 

38.5 

.0256 

-10.5 

-.1001 

-.0509 

39.5 

.0250 

-9.5 

-.1112 

-.0567 

40.5 

.0244 

-8.5 

-.1252 

-.0639 

41.5 

.0238 

-7.5 

-.1431 

-.0733 

42.5 

.0233 

-6.5 

-.1671 

-.0859 

43.5 

.0227 

-5.5 

-.2007 

-.1037 

44.5 

.0222 

-4.5 

-.2513 

-.1309 

45.5 

.0217 

-3.5 

-.3365 

-.1777 

46.5 

.0213 

-2.5 

-.5108 

-.2771 

47.5 

.0208 

-1.5 

-1.0986 

-.6479 

48.5 

.0201+ 

-.5 

0 

1.0 

49.5 

.0200 
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TABLE  II.-  COMPUTATION  B1  UNEQUAL  INTERVALS,  TRANSITION  FROM 
ONE  INTERVAL  SIZE  TO  ANOTHER 

(a)  Ax  =  0.002. 


(b)  Sx  =  0.006. 


X 

c 

% 

°n+l  - 

xn  ~  xo 

Ax 

Jno 

Jno* 

0 

.006 

.012 

.018 

.024 

a0 

°3 

a6 

a9 

CT12 

a3  - 
a6  - 
°9  ~ 
°12  " 
0i5  - 

ff0 

°3 

<J6 

a9 

a12 

-1.5 

.5 

1.5 

2.5 

-1.0986 

0 

I.O986  ! 

.5108 
.3365 

-0.6479 

1.0 

.4507 

.2338 

.1588 

.030 

.036 

.0L2 

.048 

0^4 

°15 
al  6 
a17 
sl8 

CTl6  - 
ct17  ' 

018  - 
019  - 

ct15 
°1 6 
©17 

CTl8 

3.5 

4.5 

5.5 

6.5 

7.5 

.2513 

.2007 

.1671 

.1431 

.1252 

.1204 

.0970 

.0812 

.0698 

.0613 

060 

8.5 

.1112 

.0546 

.066 

cno 

• 

• 

9.5 

10.5 

.1001 

.0910 

.0492 

.0448 

•  V  (d 

O78 

11.5 

.0834 

.0411 

.084 

12,5 

.0770 

.0380 

.090 

.096 

ff25 
ff2  6 

a26  - 
027  - 

°25 

a26 

13.5 

14.5 

.0715 

.0667 

.0353 

.0330 

TABLE  III.-  COMPUTATION  FOR  =  O.065  BY  UNEQUAL  INTERVALS 


NACA  TN  2^51 


Analytical  example  for  testing  accuracy  of  method 


\NACA 
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(t>)  Part  of  figure  4(a)  plotted  to  larger  scale. 
Figure  4.-  Concluded. 


(a)  Comparison  with  analytical  results. 

Figure  5.-  Results  obtained  by  Naiman's  method.  40-,  80-,  and  l60-point 

solutions. 


(b)  Figure  5(a)  plotted  to  a  larger  scale 
Figure  5.-  Concluded. 


(a)  Comparison  "with  analytical  results  and  results  of  Naiman’s  method 

Figure  6.-  Results  obtained  by  Multhopp's  method.  31-  and  63-point 

solutions. 
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for  example  I  by  Naiman' s  method. 
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Figure  11.-  Direct  problem  for  example  I  by  Multhopp's  method. 


__L  Solution  J _ 

Naim  an  80 -point 


78 
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Figure  12.-  Results  obtained  for  example  I  by  method  of  unequal  intervals 
compared  with  results  obtained  by  Naiman's  80-  and  160-point  solutions 
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Figure  13.-  Solution  of  inverse  problem.  Results  obtained  for  example 
by  method  of  unequal  intervals  compared  with  results  obtained  by 
Naimanfs  80-  and  l60-point  solutions. 


Solution 

o  Nctiman  80-point 
Multhopp  3J -point 


81 


Figure  15.-  Direct  problem  for  example  II  by  Multhopp's  method. 


Unequal  Intervals 

Maiman  160-point 


lire  1 6.-  Results  obtained  for  example  II  by  method  of  unequal  intervals 
compared  with  results  obtained  by  Naiman' s  80-  and  l60-point  solutions. 
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P-^e  17.-  Solution  of  inverse  problem*  Results  obtained  for  example  II 
by  method  of  unequal  intervals  compared  with  results  obtained  by 
Naiman*s  80-  and  l60-point  solutions. 


84 
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Figure  18.-  Results  obtained  for  example  III  by  method  of  unequal  inter¬ 
vals  compared  with  results  obtained  by  Naiman’s  80-  and  l60-point 
solutions. 
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